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Abstract 

^ ■ The second order perturbations in Friedmann-Robertson-Walker universe filled with 

a perfect fluid are completely formulated in the gauge invariant manner without any 
O ■ gauge fixing. All components of the Einstein equations are derived neglecting the first 

order vector and tensor modes. These equations imply that the tensor and the vector 
' mode of the second order metric perturbations may be generated by the non-linear 

effects of the Einstein equations from the first order density perturbations. 



> 



Recently, the first order approximation of the early universe from a homogeneous isotropic one is 
revealed by the observation of the CMB by Wilkinson Microwave Anisotropy Probe (WMAP)[T] and is 
suggested that fluctuations in the early universe are adiabatic and Gaussian at least in the first order 
\ approximation. One of the next theoretical tasks is to clarify the accuracy of these results, for example, 

through the non-Gaussianity. To do this, the second order cosmological perturbation theory is necessary. 

In this article, we show the gauge invariant formulation of the general relativistic second order cosmo- 
logical perturbations on the background Friedmann-Robertson-Walker (FRW) universe A4q filled with 
| the perfect fluid whose metric is given by 

^ gab = a 2 {r]) (-(di]) a (di]) b + jijidx^aidx^),,) , (1) 

5^ | where jij is the metric on maximally symmetric three space. The details of our formulation is given in 

Refs.0. 

The gauge transformation rules for the variable Q, which is expanded as Q\ = Qo + A^Q + ^X 2 ^Q, 
are given by 
rN ' 

c3 1 ^ ~ ( xQ = <2o, { yQ - ( xQ = 2£(. {1) {1 ^Q + {,££(2) + Qo, (2) 

where X and y represet two different gauge choices, and £" 2 -j are generators of the first and the 
second order gauge transformations, respectively. 

The metric g ab on the physical spacetime M.\ is expanded as g ab = g a b + Xh ab + ^-l a b- We de- 
compose the components of the first order metric perturbation h ab into the three sets of variables 
{h m , h( VL) ,h( L) , h TL }, {h(y) v h(TV)i}i and h (TT)i V which are defined by 

h vi =: Dih(y L ) + h(y)i, h {j =: a 2 h^ ij + a 2 /i (T) y, D l h {v)i = 0, 7 iJ /i( T) .,. = 0, (3) 

h(r)ij =■ ( D iDj - ^lij^j h(T L ) + 2D(ih(rv)j) + h (TT)ij D l h(TV)i = 0, D % h(TT)ij = 0- ( 4 ) 

Inspecting gauge transformation rules we define a vector field X a by 

X a := X v {dr)) a + Xi(dx % ) a , X v := h( VL) - ^a 2 d T h( TL ), , X. t := a 2 (h( T v)i + ~Dih(TL) \ , (5) 

where X a is transformed as yX a — xX a = C(i)a under the gauge transformation ([2]). We can also define 
the gauge invariant variables for the linear order metric perturbation by 

H ab = ~2a 2 ^ (di 1 ) a {di 1 ) b + 2a 2 v] {dri) (a {dx l ) b) + a 2 ^-2 * 7y+ § j (di ! )„(& 3 ') () , (6) 
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■ (!) C 1 ) ^ ■ (!) W (!) (1) (1) 

where D % Vi=X[ij]=X % i = D l Xij — 0- In the cosmological perturbations [3], {$, v%, Xij are called the 
scalar, vector, and tensor modes, respectively. In terms of the variables TL a b and X a , the original first 
order metric perturbation h ab is given by 

h a b ='■ Hab + £x9ab- (7) 

(1) (1) 

Since the scalar mode dominates in the early universe, we assume that Vi=Xij = m this article. 

As shown in Ref.[2], through the above variables X a and h a b, the second order metric perturbation 
l a b is decomposed as 

lab =■ £ab + 2£ X h a b + \£y — £ X ) 9ab- (8) 

The variables C a b and Y a are the gauge invariant and variant parts of l a b, respectively. The vector field 
Y a is transformed as yY a — xY a = + [£(1), X] a under the gauge transformations ((2]). The components 
of Lab are given by 

C ab = -2a 2 V (df]) a (df]) b + 2a 2 v} (dr])^ a (dx l ) b) + a 2 ^-2 * 7^+ Xij\ {dx l ) a {dx') b , (9) 

■ (2) (2) ^ - (2) 

where D l Ui =X[ij] =X % i= D l Xij= 0. 

As shown in Ref . \2\ , by using the above variables X a and Y a , we can find the gauge invariant variables 
for the perturbations of an arbitrary field as 

WQ:=WQ-£ x Q q ,, ^Q:=^Q-2£ X ^Q-{£ Y -£ 2 X }Q . (10) 

As the matter contents, in this article, we consider the perfect fluid whose energy-momentum tensor is 
given by T a = (e + p) u a u b + pS a b . We expand these fluid components e, p, and u a as 

(1) 1 9 (2) (1) 1 „ (2) (1) 1 (2) 

e = e + A e +-\ z e , p = p + \ p +-\ z p , Ua = u a + A u a +-Y u a p. (11) 

Following the definitions (|10p , we easily obtain the corresponding gauge invariant variables for these 
perturbations of the fluid components: 

(1) (1) „ (1) (1) „ (1) ,(*\ „ (2) (2) (1) f n n2 



£ := e -£x£, 7>:=P -j£xP, U a :=(u a ) -£ x u a , £ := e -2£ x e -{£ y -i^}e, 

(2) (2) (1) , , , (2) (2) (i) 

7> := P -2_£* P -{£ Y -£\} p , U a :={u a ) -2£ x ui - {£ Y - £ 2 X } u a . 



, , . (i) (i) 

Through g u a u b = g u a u b = —1 and neglecting the rotational part in £Y Q , the components oilA a are 

(!) (1) ' (1) 

given by U a = -a $ (dryja + oA u (dx 4 ) a ._ 

We also expand the Einstein tensor as G a b = G a b + A (1) G a b + \\ 2( - 2) G a b . From Eqs.© and ©, each 
order perturbation of the Einstein tensor is given by 

(1) G a b = (1) Q a b [H] + £xG a b , ^G a b = ^g a b [C] + ^Q b [H,H] + 2£ x ^G b + {£ Y £ 2 X } G b (12) 

as expected from Eqs. JTD]). Here, ^Q a b [H] and ^G a b [£] + ^Q a b [H,H] are gauge invariant parts of 
the frist and the second order perturbations of the Einstein tensor, respectively. On the other hand, the 
energy mometum tensor of the perfect fluid is also expanded as T a b = T a b + A' 1 ^T , a b + ^}?^ 2 'T a b and 
^T a b an d ^T a h are also given in the form 

^>T a b = ^T a b + £xT a b , ^ 2 >T a b ^^ 2 X b + 2£ x ^>T a b + {£ Y -£ 2 x }T a b (13) 

through the definitions (fT2|) of the gauge invariant variables of the fluid components. Here, ^T a b and 
( 2 ^T a b are gauge invariant part of the first and the second order perturbation of the energy momentum 
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tensor, respectively. Then, the first and the second order perturbations of the Einstein equation are 
necessarily given in term of gauge invaraint variables: 

[H] = 8nG^T a \ Wg a b [C] + WQ a h [H,H] = 8nG (14) 

_ (i) (i) 

The traceless scalar part of the spatial component of the first equation in Eq. (|14l) yields = <j> , and 
the other components of Eq. (|14p give well-known equations [3]. 

(2) _ 

Since we neglect the first order vector and tensor modes, U a is given by 

(2) //U)\ 2 (1) (1) (2)\ / (2) (2)\ 

U a = o((*J - A V D l V - $ J (drj) a + a ( A V + Vi I (dx*) a , (15) 

where D l V%= 0. All components of the second equation in Eq. (|14|) are summarized as follows: As the 
scalar parts, we have 



(2) (2) *\ / 1 

£ = (-3Hd v + A + 3K- m 2 ) $ -r + - (A~ 1 D i D j T i j - -] 

-~Hd n (A + 3K)- 1 (A 1 /)'/),! , ' ~ l^k k ^j , (16) 
87rGa 2 (e+p)A ^ = -2ft, A $ -2WA $ +AA~ 1 AT fe 

-3ft,A (A + my 1 (a '/r/;,iy - Jr fc fe ) , (17) 

(2) „.„ „ ,o\ ^2) 1 



47rGa 2 p = (ft 2 . + md v -K + 2d v U + H 2 ) $ —-A~ 1 D i DjT. 3 
3 

+ 2 

(2) (2) 



(ft 2 + 2Hd n ) (A + 3K)- 1 (A~ 1 D i D j T i j ±tA , (18) 



*-'<&' - -(A + si^r^A^AAT/'-gr^J. (19) 

where 7i = d v a/a and 

T := 87rGa 2 (e + p)L>' V A « -3A $ $ -3 15,, $ j - 8 $ A $ -12 (if + H 2 ) I $ 

, /(i) (l) (i) (i) (i) (i) (i) (i) 

:= -16nGa 2 I £ + V J A V + 12ft $ A $ -4 $ ft, A $ -4ft, $ A $, (20) 

IV := 167rGa 2 (e + p)D l ( v D j ( v -AD, $ D° $ -8 $ A A $ 

/ fc W (1) . (1) /„ (1)\ 2 , „ 9N /(1)\ 2 (1) „ U)\ , 

+2 3 A $ A $ +4 $ A $ + f ft, $ J +4 (2ft,W + if + H 2 ) ( $ J + 8ft $ ft, $ 7i J . 
As the vector parts, we have 

8nGa 2 (e+ P ) Vi = \ (A + 2JQ v\ + (r< - D i A~ 1 D k T k ) , (21) 

ft, (a 2 v^J = 2a 2 (A + 2K ) _1 {D i A~ 1 D k D l T k l - DT>} . (22) 

(2) 

As the tensor parts, we have the evolution equation of Xij 

(2) 

(ft 2 + 2Hd v +2K-A) 

= 2T« - ^r k k - 3 (DiDj - p^A^j (A + my 1 (a-^aiV - ir fc fe 

+4 (A, (A + 2K)- 1 D j) A- 1 D l D k T l k - D {i (A + 2K)' 1 D k T f)k ) , (23) 




Equations (f2"2")) and (|2^|) imply that the second order vector and tensor modes may be generated due to 
the scalar-scalar mode coupling of the first order perturbation. 

(2) 

Further, the equations (fTr?)) and (TT9"|) are reduced to the single equation for $ 

2d v H + (1 + 3c 2 s )(H 2 -K)) $ 

(2 + 3c 2 s ) Hd n ) (A + Si^)- 1 (a 1 /)'/),! , ~ ^ fc ) . (24) 
Here, we have used the second order perturbation of the equation of state for the fluid components 

(2) , (2) (2) dc 2 (1)2 „dc 2 s «« dT (I) 2 

oe ob ob 

where 5 and S are the gauge invariant entropy perturbation of the first and second order, respectively, 
we denoted that c 2 := dp/de and t := dp/dS. The equation (|24|) will be useful to discus non-linear 
effects in the CMB physics [3]. We also derive the similar equations in the case where the matter content 
of the universe is a single scalar field[2]. 

Now, we are developing our formulation to the case in which the first order vector and tensor modes 
are not negligible. In some inflationary scenario, the tensor mode are also generated by the quantum 
fluctuations. This extension is necessary to clarify the evolution of the second order perturbation in the 
existence of the first order tensor mode. Further, to apply this formulation to clarify the non-linear effects 
in CMB physics[3], we have to extend our formulation to multi-field system and the Einstein Boltzmann 
system. This extension of our formulation is one of our future works. 

Moreover, the rotational part of the fluid velocity in Eqs. ([21]) of the vector mode is also important 
in the early universe because this part of the fluid velocity is related to the generation of the magnetic 
field in the early universe[5]. The generation of the tensor mode by Eq. (1231) is also interesting, since this 
is one of the generation process of gravitational waves. We have already known that the fluctuations of 
the scalar mode exist in the early universe from the anisotropy of the CMB [6 . Hence, the generation of 
the vector mode and tensor mode due to the second order perturbation will give the lower limit of these 
modes in the early universe. 
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